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Abstract 



A set is called semidefinite representable or semidefinite programming (SDP) representable if it can 
be represented as the projection of a higher dimensional set which is represented by some Linear Matrix 
Inequality (LMI). This paper discuss the semidefinite representability conditions for convex sets of the 
form S-oif) = {x £ T> : f{x) > 0}. Here T> = {x G : gi{x) > 0, ■ ■ ■ ,gm(x) > 0} is a convex domain 
defined by some "nice" concave polynomials gi(x) (they satisfy certain concavity certificates), and f{x) is 
a polynomial or rational function. When f(x) is concave over T>, we prove that ST>{f) has some explicit 
semidefinite representations under certain conditions called preordering concavity or q-module concavity, 
which are based on the Positivstellensatz certificates for the first order concavity criteria: 



When f{x) is a polynomial or rational function having singularities on the boundary of S-oif), a per- 
spective transformation is introduced to find some explicit semidefinite representations for 5i>(/) under 
certain conditions. In the particular case n = 2, if the Laurent expansion of f{x) around one singular 
point has only two consecutive homogeneous parts, we show that 5'r>(/) always admits an explicitly 
constructible semidefinite representation. 

Key words: convex set, linear matrix inequality, perspective transformation, polynomial, Positivstellensatz, 
preordering convex/concave, q-module convex/concave, rational function, singularity, semidefinite program- 
ming, sum of squares 

1 Introduction 

Semidefinite programming (SDP) [U IIH I12[ 121] is an important convex optimization problem. It has wide 
applications in combinatorial optimization, control theory and nonconvex polynomial optimization as well 
as many other areas. There are efficient numerical algorithms and standard packages for solving semidefinite 
programming. Hence, a fundamental problem in optimization theory is what sets can be presented by 
semidefinite programming. This paper discusses this problem. 

A set S is said to be Linear Matrix Inequality (LMI) representable if 



for some symmetric matrices Ai. Here the notation X ^ ("^ 0) means X is positive semidefinite (definite). 
The above is then called an LMI representation for S. If S is representable as the projection of 



f{u) + Vf{uf{x-u)-f{x)>0, \/x,ueV. 



S = {x e K" : Ao + Aixi H h A„a;„ ^ 0} 
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that is, S = |a; G R" : 3u G R", {x,u) G S^, for some symmetric matrices Ai and then S is called 

semidefinite representable or semidefinite programming (SDP) representable. The lifted LMI above is then 
called a semidefinite representation, SDP representation or Zi/ted LMI representation for S*. Sometimes, we 
also say S = S ii S equals the projection of the lift S. 

Nesterov and Nemirovski ([11]), Ben-Tal and Nemirovski ([1]), and Nemirovsky ([H]) gave collections 
of examples of SDP representable sets. Thereby leading to the fundamental question which sets are SDP 
representable? Obviously, to be SDP representable, S must be convex and semialgebraic. Is this necessary 
condition also sufficient? What are the sufficient conditions for S to be SDP representable? Note that not 
every convex semialgebraic set is LMI representable (sec Helton and Vinnikov [Hj)- 

Prior work When 5 is a convex set of the form {x G R" : gi{x) > 0, • ■ • ,gm{x) > 0} defined by 
polynomials gi{x), there is recent work on the SDP representability of S. Parrilo [13] gave a construction 
of lifted LMIs using moments and sum of squares techniques, and proved the construction gives an SDP 
representation in the two dimensional case when the boundary of 5' is a single rational planar curve of 
genus zero. Lasserre [10] showed the construction can give arbitrarily accurate approximations to compact 
S, and the construction gives a lifted LMI for S under some algebraic properties called S-BDR or PP- 
BDR, i.e., requiring almost all positive affine polynomials on 5* have certain SOS representations with 
uniformly bounded degrees. Helton and Nie [5] proved that the convex sets of the form {x G R" : gi{x) > 
0, • • • ,gm{x) > 0} are SDP representable if every gi{x) is sos-concave {—\7^gi{x) = Gi{x)'^Gi{x) for some 
possibly nonsquare matrix polynomial Gi(x)), or every gi{x) is strictly quasi-concave on 5*, or a mixture 
of the both. Later, based on the work [5], Helton and Nie [7] proved a very general result that a compact 
convex semialgebraic set S is always SDP representable if the boundary of S is nonsingular and has positive 
curvature. This sufficient condition is not far away from being necessary: the boundary of a convex set has 
nonnegative curvature when it is nonsingular. So the only unaddressed cases for SDP representability are 
that the boundary of a convex set has zero curvature somewhere or has some singularities. 

Contributions The results in [HI [T] [TO] are more on the theoretical existence of SDP representations. 
The constructions given there might be too complicated to be useful for computational purposes. And these 
results sometimes need check conditions of Hessians of defining polynomials, which sometimes are difficult 
or inconvenient to verify in practice. However, in many applications, people often want explicit and simple 
semidefinite representations. Thus some "simple" SDP representations and conditions justifying them are 
favorable in practical applications. All these practical issues motivate this paper. Our contributions come 
in the following three aspects. 

First, there are some convex sets defined by polynomials that are not concave in the whole space R" but 
concave over a convex domain V C R". For instance, for convex set {x G R^ : X2 — x\ > {),xi > 0} , the 
defining polynomial X2 — x\ is not concave when xi < 0, but is concave over the domain R+ x R. However, 
this set allows an SDP representation, e.g., 

{xi,X2) : 3u, 

For convex sets given in the form Sv{f) = {x & V : f{x) > 0}, where f{x) is a polynomial concave over 
a convex domain V, we prove some sufficient conditions for semidefinite representability of Sxi{f) and give 
explicit SDP representations. This will be discussed in Section [2j 

Second, there are some convex sets defined by rational functions (also called rational polynomials) which 
are concave over a convex domain 2? of R". If we redefine them by using polynomials, the concavity of 
rational functions might not be preserved. For instance, the unbounded convex set 



Xi u 


h 0, 


1 Xl 




U X2 


Xi u 





x G R+ : 1 ^— > 

XiX2 



is defined by a rational function concave over R^ (R+ is the set of nonnegative real numbers). This set can 
be equivalently defined by polynomials 

{a; G R^ : X1X2 - 1 > 0, xi > 0, > O} . 
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But X1X2 — I is not concave anywhere. The prior results in [Sl[7] do not imply the SDP representability of 
this set. However, this set is SDP representable, e.g., 



X e 



Xi 1 
1 X2 



>- 



For convex sets given in the form Svif) — {x ^ T) : f{x) > 0}, where f{x) is a rational function concave 
over a convex domain J), we prove some sufficient conditions for semidefinite representability of Svif) and 
give explicit SDP representations. This will be discussed in Section [3] 

Thirds there arc some convex sets that are defined by polynomials or rational functions which arc singular 
on the boundary. For instance, the set 

{x e : xl~xl~ x\ > 0, xi > 0} 

is convex, and the origin is on the boundary. The polynomial x^—xf— x\ is singular at the origin, i.e., its 
gradient vanishes at the origin. The earlier results in [6l[7] do not imply the SDP representability of this set. 
However, this set can be equivalently defined as 

3,2 

ixi,X2) e M+ X M : xi - X-? ^ > 

Xi 

a convex set defined by a concave rational function over the domain x M. By Schur's complement, we 
know it can be represented as 





Xi 


X2 


Xi 




|(2:i,a;2) : 


X2 


Xi 





=•1 




_Xi 
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It is an LMI representation without projections. The technique of Schur's complement works only for very 
special concave rational functions, and is usually difficult to be applied for general cases. For singular 
convex sets of the form Svif) = {x € V : f{x) > 0}, where f{x) is a polynomial or rational function with 
singularities on the boundary, we give some sufficient conditions for semidefinite representability of S'-d(/) 
and give explicit SDP representations. In the particular case n = 2, we show that Sx>{f) always admits an 
explicitly constructible SDP representation when the Laurent expansion of f{x) around one singular point 
has only two consecutive homogeneous parts. This will be discussed in Section [H 

In this paper, we always assume T) ~ {x £ M" : gi{x) > 0, • • • , gm{x) > 0} is a convex domain defined by 
some nice concave polynomials gi{x). Here "nice" means that they satisfy certain concavity certificates. For 
instance, a very useful case is 2? is a polyhedra. We do not require T) or Sx>{f) to be compact, as required by 
[ni[71[TU]. When f{x) is concave over 2?, the sufficient conditions for SDP representability of S-oif) proven 
in this paper are based on some certificates for the first order concavity criteria: 

f{u) + Vf{uf{x - u) - f{x) > 0, Vx, It e P. 

Some Positivstellensatz certificates like Putinar's Positivstellensatz [16] or Schmiidgen's Positivstellensatz 
[T9] for the above can be applied to justify some explicitly constructible SDP representations for Sv{f)- 

Throughout this paper, R (resp. N) denotes the set of real numbers (resp. nonnegative integers). 
For a G N" and x £ R", denote \a\ = ai + • • • + and x" = x"^ •■•x"". B{u,r) denotes the ball 
{x £ R" : ||x — u\\2 < ?'}. A vector x > means all its entries are nonnegative. A polynomial p{x) is 
said to be a sum of squares or sos if there finitely many polynomials qi{x) such that p(x) = ^(/^(x)^. A 
matrix polynomial H{x) is called a sum of squares or sos if there is a matrix polynomial G(x) such that 
H(x) =G{xfG{x). 

2 Convex sets defined by polynomials concave over domains 

In this section, consider the convex set Sv{f) — {x £ V : /(x) > 0} defined by a polynomial f{x). Here 
T) = {x € M" : gi{x) > 0, • ■ • ,gm{x) > 0} is a convex domain. When /(x) is concave on V, it must hold 

-i?/(x, u) := f{u) + \7f{uf{x -u)- fix) > 0, Vx, w £ P. 
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The difFerence Rf{x,u) is the first order Lagrange remainder. 

2.1. q-module convexity and preordering convexity 

Now we introduce some types of definitions about convexity /concavity. Define goix) = 1. We say f{x) 
is q-module convex over V if it holds 

m I ra 

Rf{x,u) ^"^gi{x) '^gj(u)<7^j{x,u) 

i=0 \j=0 

for some sos polynomials aij{x,u). Then define f{x) to be q-module concave over V if —f(x) is q-module 
convex over V. We say f{x) is preordering convex over T) if it holds 

Rf{x,u)= (x) • • • g;;^ (x) g^^{u)---gi^]:^{u)cr^^t,{x,u) 

i^G{0,l}™ \aiG{0,1}" 

for some sos polynomials cr^,^(a:, u). Similarly, ,f{x) is called preordering concave over D if —f[x) is preorder- 
ing convex over T). Obviously, the q-module convexity implies preordering convexity, which then implies the 
convexity, but the converse might not be true. 

We remark that the defining polynomials gi{x) are not unique for the domain D. When we say f{x) is 
q-module or preordering convex/concave over T), we usually assume a certain set of defining polynomials 
gi{x) is clear in the context. 

In the special case T) = R" , the definitions of q-module convexity and preordering convexity coincide each 
other, and then arc specially called first order sos convexity. And first order sos concavity is defined in a 
similar way. Recall that a polynomial f{x) is sos-convex if its Hessian V^/(a;) is sos (see [6]). An interesting 
fact is if f{x) is sos-convex then it must also be first order sos convex. This is due to that 



fix)~fiu)^Viufix~u) 

= {x — u)'^ J f"{u 4- s{x — u)) ds dt^ {x — u) 



= {x-u)'^[^J J F{u + s{x-u))'^F{u + s{x-u))dsdtj{x~u) 
is an sos polynomial (see Lemma 3.1 of [6]). 

Example 2.1. The bivariate polynomial f{x) = x^ -\- x\x2 + xix\ -\- x\ is convex over the nonnegative 
orthant M^. It is also q-module convex with respect to M^. This is due to the identity 

i?/(x, u) = {^ui + i^i^ (^4(xi - uiY + 2{xi +X2-U1- U2Y^ + 

(^U2 + ^2:2) {^^[X2 - U2f + 2{xi +X2-U1- U2f 

2.2. SDP representations 

Throughout this subsection, we assume the polynomials f{x) and gi{x) are either all q-module concave 
or all preordering concave over T). For any h{x) from the set {/(x),5i(x), • ■ • ,gm{x)}, we thus have either 

m / m 

-Rh{x,u) = ^gj(x) Y9j{u)a^j[x,u) 
for some sos polynomials a^j{x,u), or 

i'e{o,i}'" \m6{o,i}'" 
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for some sos polynomials o-^J^^(a;, u). Let di = deg(gi)] , = \^ deg^{gi^ ■ ■ ■ g'^)'] and 
4^"^ = , ,max max deg^(gii • • • c^^,,.)! 

h£{f,gi,--- ,9m} i^e{0,l}",Aie{0,l}'" 

where deg2,(-) denotes the degree of a polynomial in x. Then define d = d)^^^^ (resp. d — dpri) when fix) 
and all gi{x) are q-module (resp. preordering) concave over T>. 
Define matrices Gq'' and c'^a^ such that 



aeN'»:|Q|<2d 
c<eN":|a|<2d 



Here mfc(2;) is the vector of all monomials with degrees < k. Let y be a vector multi-indexed by integer 
vectors in N". Then define 

N,iy)^ G^i^Vo., i = 0,l,...,m, 

0!eN":|a|<2d 

Suppose the polynomial f{x) is given in the form f{x) — X^qgN"- |Q;|<2d fa^"'- Then define vector / such 
that /^y = X^aGN": |Q|<2d fa Va- Dcfinc two sets 

i^lnodif) = {v e KaeN",|a|<2d : 2/0 = 1 , Z^?/ > 0, N,{y) h 0, VO < z < to} (2.3) 
= {y e M„gN",|a|<2d : 2/0 = 1 , /^2/ > 0, N^v) ^ 0, Vz. e {0, 1}'"} (2.4) 

via their LMI representations. Then S-dH) is contained in the image of both jO'^modif) ^^'^ ^'preif) under the 
projection map: p{y) = (2/io...o, 2/oio...o, • ■ ■ , 2/oo...oi), because for any a; e Svif) we can choose ya = x" 
such that y £ C^rnodif) and y £ C^reif)- We say S'p(/) equals C^^odif) (resp. C^reif)) if '5'i5(/) equals the 
image of C^^^^f) (resp. C^reif)) under this projection. Similarly, we say a; G C'^^odif) (^sp. x S C^reU)) 
if there exists ?/ e ^^modif) (I'csp- 2/ e Cpreif)) such that :r = p(y). 

Lemma 2.2. Assume Svif) has nonempty interior. Let {x G M" : a-^a:; = 5} be a supporting hyperplane of 
Svif) such that a'^x > b, V a; G Sxiif) and a^u ~ b for some point u G Sxiif)- 

(i) If f{x) and every gi{x) are q-module concave over T), then it holds 

m 

a^x -b - Xf{x) y^ffi(a;)g,(a;) 

i=0 

for some scalar A > and sos polynomials ai{x) such that dcgigiCi) < 2^^^^^^. 

(ii) If f{x) and every gi{x) are preordering concave over T>, then it holds 

a^x-b-\f{x)= E 9V{x)---g''^^{x)a^{x). 
!/e{o.i}™ 

for some scalar A > and sos polynomials cr^(x) such that dcg{g'^^ ■ ■ ■ g'^j^cr^ < 2dprl. 
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Proof. Since Sx>{f) has nonempty interior and the polynomials f(x),gi{x), • ■ ■ ,gm{x) are all cocnave, the 
first order optimality condition holds at u for convex optimization problem [u is a minimizer) 

min a'^x subject to f{x) > 0, gi{x) > 0, i = 1, . . . , to. 

X 

Hence there exist Lagrange multipliers A > 0, Ai > 0, Am > such that 

m 

a = AV/(u) + ^ AiVgj(u), A/(u) = Xigi{u) = • • • = Xmgm{u) = 0. 

i=l 

Thus the Lagrange function a^x — 5 — Xf{x) — X]fc=i Xigi{x) has representation 

m 

X[f{u) + \7f{uf{x ~u)- fix)') + ^ A, + \/g^{uf{x - u) - . 

i=i 

Therefore, the claims (i) and (ii) can be implied immediately from the definition of q-modulc concavity or 
preordering concavity and plugging in the value of u. □ 

Theorem 2.3. Assume V and S-oif) are both convex and have nonempty interior. 

(i) If f{x) and every gi{x) are q-module concave over D, then Svif) ~ ^'^modif)- 

(ii) If f{x) and every gi{x) are preordering concave over T>, then Svif) ~ ^^reif)- 

Proof, (i) Since Sv{f) is contained in the projection of C'^rnodU)^ o^ily need prove Sv{f) 2 ^^mod(f) ■ 
For a contradiction, suppose there exists some y G ^modif) such that x = p{ii) ^ Sv{f). By the convexity 
of Svif), it holds 

Svif)^ n |.T e M" : a^x > fe}. 

{a'^x=b} is a 
supporting hyperplane 

li X ^ S-vif), then there exists one hyperplane {a'^x = b} of S-^if) such that a^x < b. By Lemma [121 we 
have representation 

m 

a^x-b = Xf{x)+Y,9^{x)^^{x) (2.5) 

for some sos polynomials ai{x) such that deg((7,;cri) < '^d^^^^. Note that d^j^l = d. Write (T,;(a;) as 

ai{x) = md_d,(a;)^Wimd_d.(2;), i = 0, 1, . . . ,to 
for some symmetric matrices Wi ^ 0. Then the identity (j2.5p becomes (noting 



a^x-b = Xfix)+J2{9^{^>d-dMmd-dAx)'^)•W^ = Xf{x)+Y,[ J2 ^W^;" ) .W,. 

i=a i=0 \Q6N":|a|<2d 

In the above identity, if we replace each x" by , then get the contradiction 

b = f^y + Y.N,{y)»W,>0. 



T ~ 

a x 



(ii) The proof is almost the same as for (i). The only difference is that we have a new representation 



T 

a X 



b^xf{x)+ J2 ffr(a^)---C"(^K(^) 



for some sos polynomials ai,{x) such that dcg(gj'^ • • - .grn'o-jy) < 2dprl ~ 2d. Write ai{x) as 

at,{x) = md-dA^V^t^md-dAx), h 0. 
Then a similar contradiction argument can be applied prove the claim. □ 
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2.3. Some special cases 

Now we turn to some special cases about q- module or preordering convexity/concavity or SDP represen- 
tations. 

2.3.1 The q-module or preordering convexity certificate using Hessian 

The q- module or preordering convexity of f{x) over the domain 2? can be verified by solving some semidefinite 
programming. See |13[ [9] about the sos polynomials and semidefinite programming. However, in some special 
cases like T> = E" , a certificate for semidefiniteness of the Hessian V^/(x) can be applied to prove the q- 
module or preordering convexity of f{x). 

First, consider the case that gk{x) are concave over R". By concavity, it holds 

gkisx + (1 - s)u) > sgkix) + (1 - s)gk{u), V s £ [0, 1], x,ue M". 

Now we assume the following certificate for the above criteria 

gk{sx + (1 - s)u) - sgk{x) - (1 - s)gk{u) = 

(7g [X, U, S) + sal {X,U, S) + (l — SjCTj (X, U, S) + S[l — sjCTg (X, u, s) 

where ctq*''' (x, u, s), ct^'^^ (a;, u, s), 0^2'^ (a;, u, s), crg'^'' (a;, u, s) are sos polynomials in (a;, u, s). Note that the iden- 
tity (|2.6|) is always true when V is a. polyhedra, i.e., every gk{x) has degree one. 

Theorem 2.4. Suppose for every 1 < k < m, the identity \2. b]) holds. If V^/(a;) belongs to the quadratic 
module (resp. preordering) generated by polynomials gi{x),--- ,g,n{x), i.e., 

V2/(x)=^.g,(x)F,(x) \resp. V^f{x)= ^ g'^^ {x) ■ ■ ■ g'^- {x)H,{x) 
i=a \ i/efo,!}'" 

for some sos matrices Hi{x) (resp. Hu{x)), then f{x) is q-module convex (resp. preordering convex) over 
the domain V. 

Proof. First suppose V^f{x) belongs to the quadratic module generated by 51(3;), ■ • ■ ,gm{x), i.e., \7^f{x) = 
giix)Hi{x) (recall goix) ~ 1) for some sos matrices Hi{x). Then we have 

f{x)~f{u)^Vf{uf{x-u) 

= {x-u)^(^j J \/'^f{u + s{x~u))dsdt^{x-u) 



Cm „i „ 



1 \ 

gk{sx + (1 — s)u)Hi{sx + (1 — s)u) ds dt \ {x — u). 



By identity (|2.6p . we have 

nl pt 



gkisx + (1 — s)u)Hi{sx + (1 — s)u) dsdt 

Jo 
1 /-t 

/ ("'0'^'' ""i ^) + '''1'"''' ""i 'S) + (1 ~ s)a^2^ {x, u, s) + s(l — s)a'"^^ (x, u, s)jHi{sx + (1 — s)u) ds dt 
JO ^ ' 

+ 9k{x) / / sHi{sx + (1 — s)u) ds dt + gkiu) / / {I — s)Hi{sx + {1 — s)u) dsdt 
Jo Jo JO Jo 

=H^''\x,u) + gk{x)H^\x,u) + gk{u)Hi''\x,u) 

for some sos matrices h'^^ {x,u)^ {x,u), H^\x^u) (see Lemma 3.1 in So f{x) is q-module convex 
over T). 
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Second, when V^/(a;) belongs to the preordermg generated by • • • ,gm{x), i.e., 



for some sos matrices H^{x), a similar argument as above shows f{x) is preordering convex over J). □ 
Second, consider the special case that V ~ R" and the polynomial f{x) is cubic. 

" = {x : .Ti > 0, • • • ,x„ > 0} be the domain. If f{x) is a cubic polynomial concave 



Since f{x) is cubic, we have 



Theorem 2.5. Let 

over K" , then 

Proof. By Theorem 12.31 it sufBces to prove f{x) is q-module concave over R" 

-V^f{x) = Ao+ xiAi + ■■■ + x„A„ 

for some symmetric matrices Ai. When f{x) is concave over R" , — V^/(x) ^ for all x > 0. Hence all 
Ai must be positive semidefinite. This means that — V^/(x) belongs to the quadratic module generated by 



xi, ■ ■ ■ ,Xn- By Theorem 12. 4| f{x) is q-module concave over 



□ 



Example 2.6. The convex set S^2^{f ) with f{x) = 1 — {x\ + x\x2 + xix"^ + 3:2) equals 
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This is because /(x) is q-mod concave over R^ (see Example I2.ip . 

Third, consider the special case of univariate polynomials. When V = M., a. univariate polynomial is 
convex if and only if it is sos-convex, which holds if and only if it is first order sos convex. When I? = / is an 
interval, we will see that a univariate polynomial is convex over / if and only if it is q-module convex over I. 

Proposition 2.7. Let f{x) be a univariate polynomial, and I be an interval like [a,b], (—00,6] or [a, 00). 
Then f{x) is convex over I if and only if it is q-module convex over I . 

Proof. First suppose / = [a,b] is finite. f{x) is convex over [a,b] if and only if f"{x) > for all x G [a, 6], 
which is true if and only if 

fix) = (To(a;) + {x - a)ai{x) + (5 - x)c72{x) 

for some sos polynomials ctq, cri, (T2 with degrees at most 2\deg{f)/2] (see Powers and Reznick [15]). In other 
words, f{x) is convex over [a,b] if and only if its Hessian belongs to the quadratic module generated by 
polynomials x ^ a,b — x. Then the conclusion can be implied by Theorem 12.41 

The proof is similar for the case (—00, b] or [a, 00). □ 

2.3.2 Epigraph of polynomial functions 

For a given convex domain V C M", f[x) is convex over V if and only if its epigraph 

epi(/) ■.= {{x,t) eVxR: f{x)<t} 

is convex. Note that epi(/) is defined by the inequality t — f{x) > 0. If we consider t — f{x) as a polynomial 
in X with coefficients in then C^modi^ ~ /) ^^'^ ^^rei^ ~ f) ^^'^ both linear in {x, t). Therefore, if f{x) is 
q-module (resp. preordering) convex over 2?, ^^modi^^f ) (rcsp. /^^^^(t — /)) presents an SDP representation 
for epi(/). 

By Proposition 12. 7i when f{x) is a univariate polynomial convex over an interval /, we know its epigraph 
epi(/) is SDP representable and C^qj^odi^ ^ f) ^-^P representation. 
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3 Convex sets defined by rational functions 

In this section, wc discuss the SDP representation of convex set S-pif) when f{x) is a rational function while 
the domain V = {x € K" : gi{x) > 0, • • • ,gm{x) > 0} is still defined by polynomials gi. Let f{x) be a 
rational function of the form 

Here fden{x) is the denominator of f{x). We assume that /(x) is concave over the domain V. So f{x) can 
not have poles in the interior intiV) of V. Without loss of generality, assume fden{x) is positive over intiV). 
Note that f{x) is not defined on the boundary dV where fden{x) vanishes. If this happens, we think of 
Sv{f) as the closure of {x £ int{T>) : f{x) > 0}. 

3.1. The q-module or preordering convexity of rational functions 

We now introduce some types of definitions about convexity /concavity for rational functions. Let 
p{x),q{x) be two given polynomials which are positive in int{'D). We say f{x) is q-module convex over 
D with respect to {p, q) if the identity 

m / m \ 

pix)q{u) ■ Rf{x,u) = ^g,{x) ^gj{u)a,j{x,u) (3.7) 

i=0 \]=0 J 

holds for some sos polynomials <7ij{x,u). Then define f{x) to be q-module concave over 2? with respect to 
{j>,q) if —f{x) is q-module convex over T? with respect to {p,q). We say f{x) is preordering convex over T) 
with respect to (p, q) if the identity 

p{x)q(u) ■ Rf{x,u) ^ E 5r(a;)---3^-(x) [ 9'i'{^)---9^n'{^)<^^A^,^)] (3-8) 

i^eioa}" \m6{o,i}'" / 

holds for some sos polynomials a^^i^i{x,u). Similarly. f{x) is called preordering concave over T) with respect 
to (j>,q) if —fix) is preordering convex over 1) with respect to {p,q). We point out that the definition of 
q-module or preordering convexity/concavity over T) for rational functions assumes a certain set of defining 
polynomials 5i(x) for 2? is clear in the context. 

In identities p.7p or (|3.8p . there is no information on how to find polynomials p,q. However, since 
Rf{x,u) has denominator fden{x) f^i,,^{u) , a possible choice for {p,q) is 

P{x) = fden{x), q{u) = fdeniu)- (3.9) 

If the choice {p,q) in (|3.9p makes the identity (|3.7p (resp. (|3.8p ) holds, we say f{x) is q-module (resp. 
preordering) convex over V with respect to p{x), or just simply say f{x) is q-module (resp. preordering) 
convex over V if the denominator fden{x) is clear in the context. 

In the special case 2? = K", the definitions of q-module and preordering convexity over 2? coincide with 
each other, and then is called first order sos convexity when (j>,q) is given by (j3.9p . as consistent with the 
definition of first order sos convexity in Section O First order sos concavity is defined similarly. 

2 

Example 3.1. (i) The rational function ^ is convex over the domain M_|_ x M. It is also q-module convex 
over x M with respect to the denominator xi, which is due to that 

xl ul f \ , 2^2. A I , s2 

2[Xl -Ui) -\ (X2 - U2) = 2 (2;iM2 - X2U1) . 

Xi Ml V ""l ""1 / XiMJ 

4 I 2, 2 I 4 

(ii) The rational function f{x) = ^i'^^i is convex over the domain M^. It can be verified that 

/f + + + + /! + /!) + /? + + 



(/(x)-/(m)-V/(m)(x-m)) = 



(x2-t-x2)(M2-f u2)2 
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where the polynomials fi are given as below 



/i = -U1U2.T2 - uiU2xl + U1U2X2 + u\u2Xi, fe = -W2X2 + ulx2 - u\x\ + u\xi, 

/2 = -U1U2X2 + UiU2x\ + U1U2X2 ~ u\u2Xi, f-j = -2uiU2XiX2 + Uiu\xi + u\u2X2, 

/a = —^22^1X2 + u\xi — u\xiX2 + u\x2, fs = u^xl - u\x\, 

/4 = u\xiX2 — W^Xi — u\xiX2 + u\x2-, fg ~ —U1U2X1 + u\u2X2- 

/5 = ^2^2 ^ u\x2 — u\x\ + u\xi, 

So the f{x) given above is first order sos convex. □ 

Obviously, the q-module convexity implies preordering convexity, which then implies the convexity, but 
the converse might not be true. For instance, -^r^ is convex over M^, but it is neither q-module nor 
preordering convex over R\ with respect to the denominator XiX2- Note that for f{x) = it holds 

UiU2XlX2Rf(x, u) — U1U2 + X1X2U2 + X1X2U1 ~ 3xiX2UlU2. 

There are no sos polynomials (7jy_^t(x, u) such that 

ulu2XiX2Rf{x,u) = E] x'^'-X2^u'^^U2^ai,fj^{x,u). 

Otherwise, if they exist, we replace (xi, X2, ui, U2) by (x^, x^, 1, 1) and then get the dehomogenized Motzkin's 
polynomial 1 + xf 

X2 ~\~ X^X2 ~ 3x-^X2 is SOS, which is impossible (see Reznick [T7]). 

Proposition 3.2. Let CQp_q{'D) (resp. CPp^q{V)) be the set of all q-module (resp. preordering) convex 
rational functions over V with respect to {p, q) . Then they have the properties: 

(i) Both CQp^q{'D) and CPp_q{T>) are convex cones. 

(ii) If fix) e CQp^qiV) (resp. f{x) £ CPp^q{V)), then f{Az+h) e CQp,q{V) (resp. f{Az+h) G CPp,q{V)), 
where T) = {z : Az + b € P} and p{z) = p{Az + b), q{z) ~ q{Az + b). That is, the q-module convexity 
or preordering convexity is preserved under linear transformations. 

Proof. The item (i) can be verified explicitly, and item (ii) can be obtained by substituting Az + b for x and 
noting the chain rule of derivatives. □ 

3.2. SDP representations 

Now we turn to the construction of SDP representations for Svif). Recall go{x) = 1. Throughout this 
subsection, we assume the polynomials f{x) and gi{x) are either all q-module concave or all preordering 
concave over V with respect to {p,q). Thus for any h{x) from {f{x),gi{x), ■ ■ ■ ,g,„(a;)}, we have either 

m / m * 

-p{x)q{u)Rh{x,u) = ^gi{x) ^ gj{u)a'lj{x,u) 
for some sos polynomials cr^j{x, u), or 

-p{x)q{u)Rh{x,u)= E 9? {x) ■ ■ ■ g'^n {x) \ E 3i'(")-".9m"Hc^!^,p(a;,w) 
for some sos polynomials ^{x,u). Let 



dirl= max max dcg^Cffi' ' ' ' 5m"cr^,^.)l ■ 

h€{f,gi,--- ,grr,} I'eio.il'^.Aieio,!}'" 
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Then set d = d^J^^^^ (resp. d = dp^l) when f{x) and gi{x) are all q-module (resp. preordermg) concave over 
V with respect to {p, q). 



Define matrices Pa'^ , Pi'^'' , ' , Q^^' such that 



Pixf ^d-d^ {x)md-d^ {xY 



< i < m, 



ieN":|a| + |L£;(p)|<2rf peN":P<LE{p) 



)j£_ 

p{x) 



Q6N":|a| + |L_B(p)|<2d ,36N":/3<L£;(p) 

(3.11) 

Here LE{p) denotes the exponent of the leading monomial of p{x) under the lexicographical ordering {xi > 
X2 > ■ ■ ■ > Xn)- Note that the miion 

{x" :aeN",|a| + |L£;(p)| < 2d} U : /3 £ N", /3 < 

is a set of polynomials and rational functions that are linearly independent. 

Let y be a vector indexed by a e N" such that |q;| + \LE{p)\ < 2d, and z be a vector indexed by /3 G N" 
such that P < LE{p). Then define 



aeN":|Q| + |L£;(p)|<2d 
aeN":\a\ + \LE{p)\<2d 



Qo^Ua 



Qa ya 



E P^'^zp^ 0<i<m, 

PeN":f3<LE(p) 

h E {0,1}" 

/3eN":/3<L£;(p) 



(3.12) 



Suppose the rational function f{x) is given in the form 



QeN":|QH-|L£;(p)|<2d 

then define vectors f^^\ Z*^^' such that 



l3m^:l3<LE{p) 



QeN":|Q| + |L_B(p)|<2d /3GN":/3<L_B(p) 



Define two SDP representable sets 



"/(I)- 


T 


2/ 


/(2)_ 






"/(I)- 


T 


y 


_/(2) 







> 0, Q^iy,z) ^ 0, VO < i < m 

>o, QAy,z)hO,yiye{o,iy 



(3.13) 
(3.14) 



We say S-pif) equals Ti^rnodif) (resp. 7^^re(/)) if '5'p(/) equals the image of Tl^^^dU) (resp. Tl^reU)) under 
the projection p(y, z) = (yio...o, yoio...o, • ■ • , yoo...oi)- Similarly, wc say x G 7^^„^^(/) (resp. a; £ Tl^reif)) 
if there exists (y, z) G Tl^rnodif) (^esp. (y, 2) G Upreif)) such that a; = p(2/,z). 

Lemma 3.3. Assume D and Svif) are both convex and have nonempty interior. Let {x G M" : a^x = b} be 
a supporting hyperplane of S-pif) such that a^x > b for all x G Sx>{f) and a^u = b for some point u G Sx>{f) 
such that q{u) > 0, and either f{u) > or fden{u) > 0. 

(i) If f{x) and every gi{x) are q-module concave over V with respect to {p,q), then 

m 

p{x) ■ (a^x -b- Xf{x)) ^^gi{x)ai{x) 

1=0 



for some scalar A > and sos polynomials ai{x) such that deg{giai) < 2d^I2.od 
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(ii) If f{x) and every gi{x) are preordering concave over V with respect to {p,q), then 

p{x)-{a^x-h-Xf{x))^ J2 g'(H^)---9!^{^M^) 

for some scalar A > and sos polynomials cri{x) such that dcg(.gj^ • • ■ g^^cr^) 

Proof. Since Sxi{f) has nonempty interior, the first order optimality condition holds at u for convex opti- 
mization problem (u is one minimizer) 



T 

mm a x 



subject to f{x) > 0, gi{x) > 0, i = 1, 



If f{u) > 0, the constraint f{x) > is inactive. If fden{u) > 0, f{x) is differential at u. Hence, in either 
case, there exist Lagrange multipliers A > 0, Ai > 0, . . . , A,„ > such that 

m 

a = AV/(ii) + ^ XiVgi{u), Xf{u) = Xigi{u) = ■■■ = A,„.g„i(ii) = 0. 

i=l 

Hence we get the identity 

m 

a^x - b- Xf{x) - ^ A,5j(x-) = 
1=1 

m 

a(/(u) + Vf{uf{x -u)- fix)) + X,(g,{u) + Vg^{uf{x ~ u) - g,{x)) . 

i=l 

Therefore, the claims (i) and (ii) can be implied immediately by the definition of q-module concavity or 
preordering concavity of / and gi, and plugging the value of m. □ 

Theorem 3.4. Assume V and S-pif) are both convex and have nonempty interior. Let {p,q) be given in 
{531 or fSJj. Suppose dim(Z(/) n Z{fden) H dSv{f)) < n - 1 and dimiZ{q) D dSvif)) <n-l. 

(i) If f{x) and every gi{x) are q-module concave over T> with respect to (p, q), then Svif) ~ T^modif) ■ 

(ii) // f{x) and every gi{x) are preordering concave over T> with respect to (p, q), then Svif) = ^^re(/)- 

Proof, (i) Since S-vif) is contained in the projection of 7^^mod(/)j we only need prove S-pif) ^ T^modif)- 
For a contradiction, suppose there exists some (y, z) g T^modif) ^VlC\\ that x = p{y) ^ S-oif). By the 
convexity of Sxiif), it holds 

Sv{f)= n {x e M" : a^x > &}. 

{a^x=6} is a 
supporting hyperplane 

If i; ^ Svif), then there exists one supporting hyperplane {a-^x = b} of S-pif ) with tangent point u G dSvif) 
such that a-^i < b. Since dim(Z(/) n Z{p) n 9S'x)(/)) < ?i — 1 and dim(Z(g) n dSx>) < n — 1, by continuity, 
we can choose {a^x = b} such that cither f{u) > or p{u) > 0, and q{u) > 0. By Lemma [3.31 we have 

a^x-6 = A/(x)+V^a.(x) (3.15) 

for some sos polynomials cri{x) such that deg{giai) < 2^^^''^^. Note that d = d^^^^. Then write ai{x) as 

ai{x) = md^d,{x)'^Wim.d-di{x), i = 0, 1, . . . ,to 
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-0.5 







0.5 



Figure 1: The convex set defined by xf + X2 > xf + XIX2 + X2 



for some symmetric matrices Wi >z 0. Then identity (|3.15p becomes (noting (|3.1ip ) 



T 

a X 



b = Xf{x) + J2 



1=0 



p{x) 



■md~d,i'^)^d-diix) ]»Wi 



i=0 \a<£N^:\a\ + \LE{p)\<2d 



peN":l3<LE{p) 



p{x) 



•W,. 



In the above identity, if we replace each x*^ by y^i -^-^ by z^, then we get the contradiction 

i=0 

(ii) The proof is almost the same as for (i). The only difference is that 



T 

a a; — 7 



A/(a;) 



E 

i/e{o,i}' 



9? 



p{x) 



-a^{x) 



for some sos polynomials cr^(x) such that dcg{giai) < 2dp^l. Note that d ~ dp^l. 



A similar contradiction 
□ 



argument like in (i) can be applied to prove the claim. 

Example 3.5. The convex set {x G M.'^ : x^ + x\ > x\ + xfx"^ + X2} can be defined as 8^2 (/) with rational 
4 1 2 2 I 4 

function /(x) = 1 - ^ . The set S'k^C/) is the shaded area bounded by a thick curve in Figure [T] 

We have already seen that f{x) is first order sos concave. So 8^2 {f) = T^jnodif)- ^ polynomial division 
shows that ^a^^a VCI2 {x)m2 {x)'^ equals 












1 





0" 




1 


Xl 


X2 


-xl 


XIX2 


xl 





1 





Xl 


X2 







Xx 


-xl 


X\X2 


-xixl 


-xl 


xixl 











X2 








^ 1 


X2 


X\X2 


X2 


-xl 


xixl 


xl 


1 


x-i 


X2 


2 2 
Xl - X2 


XIX2 


X2 


p{x) 


—X2 


—x\x\ 


— X2 


xi 


-xixl 


-xi 





X2 





X1X2 


xl 







X1X2 


-xi 


xixl 


-X1X2 


-xi 


X1X2 











xl 





0_ 




. 4 


2 

2:1X2 


xl 


-xt 


xixl 


xi 
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So we can see that S'r2(/) = 7^^„jojj(/) can be represented as 



' ( 


Xi 


X2) 


G M 


X M : 3 


2/jj 




s.t 


1 > 


2/20 ^ 


- 204, 












"0 








1 





" 




^00 


zio 


zoi 


-Zq2 


^11 


202 









1 





Xl 


a:;2 







zw 


-Z02 


Zll 


~Zi2 


-203 


212 















X2 








+ 


zoi 


Zll 


Z02 


-203 


Zl2 


203 






1 


Xi 


X2 


2/20 " 2/02 


2/11 


2/02 




-Z02 ' 


-Z12 


-Z03 


zoi 


-213 


-204 














2/11 


2/02 







Zll 


-Z03 


Z12 


-Zl-i 


-204 


213 















2/02 





_ 




_ Zq2 


Z12 


Z03 


-Z04 


213 


204 _ 


> 



The plot of the projection of the above coincides with the shaded area in Figured] 



3.3. Some special cases 



3.3.1 Epigraph of rational functions 

The rational function j[x) is convex over the convex domain 2? if and only if its epigraph 

epi(/) := {(x, e P X K : /(x) < /} 

is convex. The LMI Ti^^^^it — /) and VJ^^J^t — /) can be constructed by thinking of t — /(x) as a polynomial 
in X with coefficients in t. So, if j{x) is q-module (resp. preordering) convex over P, Tl^modi^ ^ /) (resp. 
7?.^g(i — /)) is an SDP representation for epi(/). 

Now we consider the special case that f{x) = is a univariate rational function convex over an interval 
/ = [a, b] and p{x) is positive over (a, b). Note that I = {xGM.: x~a>0,b — x>0} (x — a>Ois not 
required if a = —00, and similarly for b — x > 0). 

Theorem 3.6. Let f{x) ~ be a univariate rational function and p{x) is a polynomial nonnegative over 
an interval I. If f{x) is convex over I, then its epigraph epi(/) = T^lmodif)- 

Proof. First assume / = [a, b] is finite. Since T^^q„ilaif) is linear in t, it sufRces to show for any fixed t 

{x e [a, 6] : 3y, z, s.t., x = yi, {y, z, t) e T^^grnldif)} = ^ [«' ^1 ■ fi^) ^ 0- 

In the above the right hand side is contained in the left hand side. Now we prove the converse. For a 
contradiction, suppose there exists a tuple (y, 2,t) e T^qmod^f) s^ndti that and x = yi does not belong to the 
convex set {x G [a, b] : f{x) < t}. Then there exists an affine function cq + cix such that 

Co + cix < 0, Co + cix > 0, y X G [a,b] : f{x) < t. 

Since f{x) is convex over [a,b], by optimality condition, there exists A > such that 

„/ ^^ (Cn + C-\X — Xt)p(x) + Xq(x) r .-, 

co + cix^X{t-f{x))=^-^ ^ ^>0, yxe[a,b]. 

p{x) 

Then we can sec that (co + cix — Xt)p(x) + Xq{x) is a univariate polynomial nonnegative on [a, 6]. So 
there exist sos polynomials sa{x), si{x), S2{x) of degrees at most 2d, 2d — 2, 2d — 2 respectively (assume 
deg(p) + 1, dcg{q) < 2d and see Powers and Reznick [15]) such that 

(co + Cix — Xt)p{x) + Xq{x) = Sq{x) + (x — a)si{x) + (6 — a)s2{x) 

and hence then 

CO + cix = X{t - fix)) + ^oi-) + i^-^)^ii^) + ib-a)s2{x)_ 
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Write so{x), si{x), S2{x) as 

so{x) = md{x)'^Womd{x), si{x) = md-i{x)'^Wimd-i{x), S2{x) = md~i{x)'^W2m.d-i{x) 
for some symmetric Wq, Wi, W2 ^ 0. Then we have (noting go ~ 1, gi ~ x ~ a, g2 = b — x and p. lip ) 

ax *• — I X I 



!=0 

2 



-b = Xit- fix)) + ( ^md-d^md-d. (x)^ ) • (do = 0, di = = 1) 



= x{t-f{x)) + Y^{ g«.'^+ 5: ^'^^l'^- 

'i=0 yaeN:|Q| + |L£;(p)|<2d /3GN:/3<L£;(p) 'J 

In the above identity, if we replace each x" by Uq, and -7^ by zg, then get the contradiction 

2 

a^x - 6 = A(t - (/(i))^y - (/(2')^z) + Q^iv^ ^)'W,> 0. 

1=0 

Therefore we get epi(/) = 7e[°;^^(/). 

When / is an infinite interval of the form [a, 00) ~ {x : x — a > 0}, (—00, 6] = {a; : & — a; > 0} or 
(— cx), Qo) = {a; : 1 > 0}, a similar argument can be applied to prove epi(/) = Ti-qmodif)- ^ 

3.3.2 Convex sets defined by structured rational functions 

For the convenience of discussion, we define some basic convex sets (r < s) 

Kr^s = {{w, v)eR\xM.+ : wi---Wr> v'} (3.16) 

which are all SDP representable (see §3.3 in P). 

First, consider epigraphs of rational functions of the form 

f{x) = (3.17) 

Pl{x) ■ --Prix) 

where qi{x) , pj (x) are polynomials nonnegative over T) and the integer Si > r + 1. Then 

epi(/) - |(x,i) ePxK: Yiq,{x))'' < pi(x) • ■ •p,(x)t| . 
The epigraph epi(/) can be equivalently presented as 

epi(/) = < (a;, i) e P X M : 3u„ Wi, Y — — ^ *(a;) < Ui, m < Piix), i = 1, . . . , P > . 

Note that 

q,{x) < Ui (x, Mj) e epi(gi), < pi{x) (x, ~Wi) G epi(-p,), 

P ^s, 

^ < t • • ■ , Wp > 0, Ui H hVp <t,{wi,--- ,Wr,Vi,Ui) e Kr+l,s,- 

'-^ W\ - ■ -Wr 

Hence we obtain that 

p 

epi(/) = |(a;,t) S I? x K : ~3u^,w^,Vi, ^ < t, (wi, • • • ,Wr,Vi,uC) G Kr+i.si, 

i=l 

{x,Ui) g epi(gj), (a:, -Wi) £ epi{-pi),i = 1, . . . 

When gi(x) are q-module (resp. preordering) convex over V, we know epi{qi) = C^j^odi^i " (resp. 
epi(gi) = C^^^{ui — qi)), and similarly for —pi{x). Thus, we get the theorem: 
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Theorem 3.7. Suppose f{x) is given in the form l[3.17\ l and all qi{x),pj{x) there are nonnegative over T). 
(i) If all qi{x), —pj{x) are q-module convex overT), then 

p P 



epi(/) = [{x,t) : 3ui,Wi,Vi,x e Q Cf„^odiui - Qi) f] ^fmodiPj - Wj), 

i=l 3=1 

P 

< t, (Wi, • ■ • ,Wr,Vi,Ui) S Kr+l.s}^. 

i=l 

(ii) // all Qi (x) , —pj (x) are preordering convex over T), then 

p P 

epi(/) = \{x,t) : 3u,,Wi,V,,X £ f] ^prd'^i ~ Qi) f] ^preiPj - Wj), 

1=1 j=i 

p 

y^^Vj < t, {Wi, ■ ■ ■ ,Wr,Vi,Ui) £ Kr+l,s,y 



Example 3.8. (i) Consider the epigraph ^^{xi,X2,t) €R+ xRxR:t> -^^^^^j 
here is given in the form (|3.17p . By Theorem 13.71 its epigraph can be represented as 



The rational function 



{xi,X2,t) : 3u, 



t u 

U Xl 



^0, 



U — 1 X2 
X2 1 



>- 



(ii) Consider the epigraph G i?(0, 1) x R : t> (il^J)^.^(\_^2) | • The rational function here is given 

the form (|3.17p . By Theorem 13.71 its epigraph can be represented as 



(a;, i) : 3 u, w^, wi, s.t., ^ < t, 



i=l 



1 — Wi Xi 
Xi 1 



^0, 



U X 
X In 



^ 0, (u, U)i, . . . , W„, u) e Kn+ 



l.n+1 



Second, consider epigraphs of rational functions of the form 

L 



(3.18) 



fc=i 



where fk{x) are rational functions given in form (j3.17p and > 1. Then the epigraph epi(/i) = {{x,t) G 
2? X R : h{x) < t} can be presented as 

epi{h) = |(x,i) : 3rik,Tk, ^ < {x,r]k) G epi(/fc), {Tk,r]k) G I<i^b^, fc = 1, • • • ,l| . 

Once the SDP representation for each epi(/fc) is available, one SDP representation for epi(/i) can be obtained 
consequently from the above. 



sion, we know it can be represented as 

{xi,X2,t) : 3u, 



t u 
u 1 



^0, 



R+ 


X R : 


t > 


( 1+^2 


Xl 





r 







Xl 


X2 


"1 


1 


X2 


U 
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Third, consider the convex sets given in the form 

T = jx e P : a{x) > ^ hk{x) 



(3.19) 



fe=i 



where a(x) is a polynomial and every hk{x) is given in the form (|3.18p . Then 



T = 



ja; : 6*1, ■ • • , 6*7, t > ^ ffc, (a;, -t) e epi(-a(a;)), {x, Ok) £ eTpi{hk{x)), = 1, . . . , j| 



When a{x) is q-modulc or preordcring concave over D, epi(— a(x)) is representable by '^^^^(0(2;) — t) or 
C^^^{a{x) — t). Once the SDP representations for epi(— 0(0;)) and all epi(/),fc(x)) are available, an SDP 
representation for T can be obtained consequently. 



Example 3.10. (i) Consider the convex set T = |(a;i,a;2) G 



1 - ixl > 



(1+ 



the form ((3TT9)) with V 



T 



{xi,X2) : 3u,w, 



From the above discussion, we have 

- 1 X2 



ill} 



It is given in 



1 — U X2 
X2 i 



^0, 



u w 

W Xl 



X2 1 

>- 0. 



l-xl> 



(1+ 



Note that [u^xi,w) G -f5r2,2 has the representation 

(ii) Consider the convex set T = |(xi,X2) G R+ x 
I> = K.+ X R. From the above discussion, we know it can be represented as 

[1 + ^ 



1 2.2-12 ~| 

./ > . It is given in the form (I3.19P with 
^1 J 



(a;i,2;2) 



X2 



X2 

1 - u 



^0, 



Xl 





1 ' 







Xl 


X2 


=•1 


1 


X2 


U 





4 Convex sets with singularities 

Let Sv{f) = {x <E V : f{x) > 0} be a convex set defined by a polynomial or rational function f{x). Here 
V = {x E M" : gi{x) > 0, • ■ • ,gm{x) > 0} is still a convex domain defined by polynomials. Suppose the 
origin belongs to S'p(/) and is a singular point of the hypcrsurface Z{f) = {a; G C" : f{x) ~ 0}, i.e., 

/(0,0) = 0, V/(0,0)=0. 

We are interested in finding SDP representability conditions for Sv{f)- 

As we have seen in Introduction, one natural approach to getting an SDP representation for St> (/) is to 
find a "nicer" defining function (possibly a concave rational function). Let p{x) be a polynomial or rational 
function positive in int{T)). Then we can see 5p(/) is the closure of the set 

X G int{V) : > 

p{x) 

If has nice properties, e.g., has special structures discussed in Section |3l or it is q- module or 
preordering concave over T), then an explicit SDP representation for S-pif) can be obtained. For instance, 
consider the convex set 

{{xi,X2) G R+ : -xl + Sxixl - (a-^ + xlf > 0}. 
The origin is a singular point on its boundary. If we choose p{x) = X2, then it can be presented as 

„3 ^4 • 



f a;^ a;^ 1 

<^ (xi,a:2) gR^ : 3xi>2xl+xl + -^ + -^\ . 
I 2;2 a;2 J 
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Then this set can be represented as 

|(a;i, a;2) G : 3 ui, W2, U3, 1*4, s.t., 3a;i > 2ui + U2 + U3 + "4, (ui, Xi) G i^i,2, 
iu2,X2) e A'1,2, (w3,a;2,a;2,a;i) G ^3,3, (u4, a;2, z^2, a;i) G A'3,4|. 

However, there is no general procedure to find such a nice function p{x). In convex analysis, there is a 
technique called perspective transformation which might be very useful now. Generally, we can assume 

Sv{f)cR+xW'-\ int{Sv{f))^9. 

Define the perspective transformation p as 

p(a;i,a;2, • • • ,a;„) = (l/xi, a:2/xi, • ■ • ,Xn/xi). 

The image of S-pif) under the perspective transformation p is 

{{1/XUX2/XU--- ,Xn/xi) : xeSvif)} CR+ xR"-\ 
which is also convex (see §2.3 in [5]). Define new coordinates 

^ _ Xn 

X\ — , X2 , * ' * 7 Xji — . 

Xl Xl Xl 

Denote x = {x2, • ■ ■ Suppose f{x) has Laurent expansion around the origin 

f{x) = fb{x) - fb+i{x) fd{x) 

where every fk{x) is a homogeneous part of degree k. Let 

fd-b{x) 



m ■■= Mi) 



where fi{x2, ■ ■ ■ := x\'^^ fb-i-i{x) . Define a new domain 2? as 

P = {iGM" : >0,--- ,g,n{i)>0} 

where gi{x) = -as^^-y. Note that V is convex if and only if V is convex (see §2.3 in [5]). Therefore, under 

^1 

the perspective transformation p, the set Sv{f) can be equivalently defined as 
Proposition 4.1. If Sv{f) is convex, then fo{x) > for any x — {xi,x) G Sff{f). 

Proof. Fix X ~ {xi, • ■ • , Xn) G Sv{f) and x = (xi, • • ■ , x„) G Sfi{f) such that x — p(xi, X2, ■ • • , a;„). By the 
convexity of Sv{f), the line segment {tx ; < t < 1} belongs to Sv{f)- Thus its image 

lp{tXi,tX2, ■ ■ ■ ,tXn) = {^Xi,X2, • ' • , Xn) 

belongs to Sf,{f). Now let i ^ 0. Then f{jXi,X2, ■ ■ ■ ,in) > imphes /o(a;) > 0. □ 

4.1. The case of structured fk{x) 

In this subsection, we assume fk{x) have special structures. Then the methods in Subsection 13.3.21 can 
be apphed to construct SDP representations for Sf,{f). 
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Theorem 4.2. Suppose every fk{x) (fc = 1, ■ • • ,d — h) is given in the form 

i=i j=i i=i 

for some polynomials qk^i{x),pk,j.e{x) which are nonnegative over T) and integers rk.i > k + I, Sk.jj > 2 
(pk,j{x) can be any affine polynomial when Skj/ is even). Then Sf,{f) can be represented as 



d—b I Lfc 



X : 3u,Uks,Vk,],£, > X! 



{uk/r 



k=i \ 1=1 



V 



Xi 



) 



e epi(-/o), {x,Uk,t) £ epi{qk,r), {x,Vk,j,i) G epi(pfej,£)|. 

Furthermore, if —foix) o-nd all qk.j{i)TPk,j.i{x) are q-module or preordering convex overT), then Sf,{f) is 
SDP representable. 

Proof. The first conclusion is obvious by introducing new variables u,Uk.i,Vk.j,e- Note that 



d—b Lk I \r-, ■ Rk Qkj , \<^i n I. 

fc=i 4=1 1 j=i e=i 



Xi 



is equivalent to 



d-b Lfc 



fc=l i=l 



(^fe,Cfc) e Ki^k, {Ck,Xi,Vk,j,e) G K2,sk,j,e- 



When — /o(i) is q-module (resp. preordering) convex over V, {x, —u) G epi(— /q) is representable by 
^%nod{'^ ~ fo) (rcsp. C'^j.^{u — /o))- Similar results hold for qkj{x),pk,j,e{x). Once the SDP representa- 
tions for epigraphs of —fo{x), qk.j{x) and pkjj{x) are all available, we can get an SDP representation for 
Sf,{f) consequently. □ 

Now we show some examples on how to apply Theorem 14.21 



Example 4.3. (a) Consider convex set S = {(xi, X2) G M+ x M : x\ — x\ — x^> 0}. Its boundary is a cubic 
curve and the origin is a singular node. This convex set is the shaded area bounded by a thick curve in 
Figure [2^a). The thin curves are other branches of this cubic curve. After the perspective transformation, 
we get 

(ii,i2) gM+xM: 1-^2-— >0 



Xi 



which can be represented as 



(ii, X2) '. 3 M, s.t.. 



\ — u in 

X2 1 



^0, 



u 1 
1 Xl 



y 



After the inverse perspective transformation, we get 



(xi,X2) : 3w, s.t.. 



Xl — U X2 
X2 Xil 



U Xl 
Xl Xl 



y 



The plot of the projection of the above coincides with the shaded area in Figure [^IJa). 

(b) Consider convex set S = {{xi, X2) G M+ x R : x\— x\— X2>(^}. The origin is a singular tacnode on the 



19 



-0.6 -0.4 -0.2 0.2 0.4 0.6 



(a): xl - xl - xl > 0, xi > 



-0.6 -0.4 -0.2 0.2 0.4 0.6 



(c): xl - Sxixj ~ {xl + x\f > 0, xi > 



-0.4 
-0.6 

-0 



(b): x\ - x\ - x\>^,Xi>Q 



C •• 



(d): xl-xl~ {xl + xlf > 0, Xl > 



Figure 2: Four singular convex sets discussed in Example 14.31 



boundary. This convex set is the shaded area bounded by a thick curve in Figure [2l^b). The thin curve is 
the other branch of the singular curve x\ ~ X2 — x\ = 0. After the perspective transformation, we get 



S = <j {XI,X2) G M+ X 

which can be represented as 

{xi,X2) ■ 3ui,U2, S.t 

After the inverse perspective transformation, we get 
b = { (Xi,X2) : dwi,M2, s.t., 

U2 Xl + Ul 



1 > 



1 + 



1 - Ul U2 
U2 l+Ui 



h 0, 



Ul 1 
1 Xl 



h 0, 



U2 X2 
X2 Xl 



>- 



^0, 



Ul Xl 
Xl 1 



^0, 



U2 X2 
X2 1 



>- 



The plot of the projection of the above coincides with the shaded area in Figure [^b). 

(c) Consider convex set S = {{xi,X2) £ K+ x R : xf — Sxix^ — {xf + x^)^ > 0}. The origin is a singular 
point on the boundary. This convex set is the shaded area bounded by a thick curve in Figure [H^c) . The 
thin curves are other branches of the singular curve xf — Sxixf — {xf + x^)"^ = 0. After the perspective 
transformation, we get 



1-3x1- 



> 



Xl 
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which can be represented as 



(ii, 2:2) : 3u,w, s.t., 



1 — U X2 
X2 i 



h 0, 



u w 

W Xi 



h 0, 



lU — 1 £2 
X2 1 



^ 



After the inverse perspective transformation, we get 



S* = < {xi,X2) '■ 3u,w, s.t. 



Xi — U X2 

1 



X2 



Xi 



u w 
w 1 



^0, 



W — Xi X2 
X2 Xi 



>- 



The plot of the projection of the above coincides with the shaded area in Figure [^c). 



(d) Consider convex set S = {{xi,X2) S IR+ 



^2 — [xf + 2:2)^ > 0}. The origin is a singular point 



on the boundary which is one branch of the lemniscate curve xf — X2 — {x\ + x^Y = 0. This convex set is 
the shaded area bounded a thick lemniscate curve in Figure [2jd). The thin curve is the other branch of the 
lemniscate curve. After the perspective transformation, we get 



{xi,X2) e M+ X M : {l- xl)> 



Xl 



which can be represented as 

{xi,X2) ■ 3u,w, s.t. 



1 + U X2 
X2 1 — U 





u 


1 


X2 






1 


Xl 









_X2 





Xl_ 





After the inverse perspective transformation, we get 



Xi+ u 
X2 



X2 



(xi,X2) : 3w, s.t. 

The plot of the projection of the above coincides with the shaded area in Figure [^d). 







u 


Xl 


X2 








Xl 


1 





=•1 


u 





1 






X2 





4.2. The case of two consecutive homogeneous parts 

In this subsection, we consider the special case that f{x) ~ fb(x) — fh+i{x) having two consecutive 
homogeneous parts. Then, after perspective transformation, we get 

5p(/) = |ie^:/o(5)-^>o|. 

By Proposition 14. 11 for any x E V, we have fo{x) > 0. Let V be the intersection of {x : fo{x) > 0} and the 
projection oiV CR+ x M"-i into M""!. Then we get 

fo{x) 

Then Sf,{f) is the epigraph epi(/i) of the rational function h{x) over V . Note that h{x) is convex over the 
domain 2?' if and only if 5*^5 (/) is convex, which then holds if and only if S-pif) is convex. So, if h{x) is 
q-module or preordering convex over V , then Ti^jnodi^i ~ ^) T^preixi — h) is an SDP representation for 
epi(/i), and then one can be obtained for Sx){h) after the inverse perspective transformation. 
A very interesting case is n = 2. Then D' must be an interval / of the real line. 

Theorem 4.4. Let n = 2 and V = I be an interval as above. If Sx>{f) is convex, then Sf,{f) = Ti-qmodi^i ~ 
h), and hence Sv{f) = P"H^gmod(*i " ^))- 
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Proof. When n = 2, f{x) is a univariate rational function. When Sv{,f) is convex, Sf,{f) is also convex. 
Since Sf,{f) is the epigraph of /(i), /(i) must be a univariate rational function convex over the interval /. 
By Theorem 13.61 its epigraph is representable by 'R-q^nodi^^ ~ Thus Sf,{f) — T^l„iod(^i ~ After the 
inverse perspective transformation, we can get an SDP representation for Sx){f)- d 

Now we see some examples on how to find SDP representations for singular convex sets by applying 
Theorem 14.41 

Example 4.5. (i) We revisit the singular convex set (a) in Example 14.31 After the perspective transforma- 
tion, we get S = |(ii,i2) G M x [—1, 1]: xi> ^^-^ | , which can be represented as 





3 2/2,2/3,^4,^0 




Xl 


> Zq 




■ 1 

X2 


- 2/2 X2 

- 2/3 2/2 


- 2/3 

- 2/4. 




< ixi,X2) 




1 X2 2/2 




zo 




Z\ 




-1 + 


zq 


> 






X2 2/2 2/3 


h 0, 


Z\ 




-1 + 


zq 


-~X2 4 


Z\ 


h 






_2/2 2/3 2/4_ 




-1 + 


Z{) - 


-X2 + 


Z\ 


-2/2 - 1 


+ 





Applying the inverse perspective transformation, we get an SDP representation for S 



3 2/2,2/3,2/4,^0 


^1, 




1 > ^0, 




X\ 
X2 


- 2/2 a;2 

- 2/3 2/2 


- 2/3 

~ y^. 


^0, 


< {xi,X2) ■■ 


Xl X2 2/2 




Zo 








-Xl 4 


Zo 


> 




X2 2/2 2/3 




Z\ 


-Xl 4 


Zo 


-X2 4 


Z\ 


> 




.2/2 2/3 2/4. 




~x\ 4 


- Zo -a;2 4 


Z\ 


-2/2 - Xl 


4zo_ 


> 



Interestingly but not surprisingly, the plot of the above coincides with the shaded area in Figure [^l^a). 
(ii) Revisit the singular convex set (c) in Example 14.31 After the perspective transformation, we get S 



|(ii,X2) G 



^[-l,l]:5i>-f - 
3 2/2,2/3,2/4,^0,2:1, 



- 16 
9(1-3^;^) 

X\ > 



I , which equals 



\V2 



16 , 



(xi,i2) : 



" 1 


i2 


2/2 




X2 


2/2 


2/3 




.2/2 


2/3 


2/4. 





Zo zi 
zi i(zo-l) 

i(zo - 1) h{zi~X2) 



1 - 32/2 X2 - 32/3 
X2 - 32/3 2/2 - 32/4 
i(^o-l) 
i(zi - i2) 
^2/2 4 |(zo - 1) 



3 

Applying the inverse perspective transformation, we get an SDP representation for S 



3 2/2, 2/3, 2/4, 2^0, 2;i, 



1 > -h2 



Ixi 



-qZo, 



ixi,X2) 



Xl 


X2 


2/2 




X2 


2/2 


2/3 


>= 0, 


.2/2 


2/3 


2/4. 





Xl - 32/2 X2 - 32/3 
a;2 - 32/3 2/2 - 3274 
zi i(zo-xi) 
|(zo-a::i) 3(^1-2:2) 
i(zo-xi) 3(21 -X2) -^2/2 4 |(zo - Xl) 



Zo 
zi 



h 
y 



> . 



Also interestingly but not surprisingly, the plot of the above also coincides with the shaded area in Figuredfc). 

Now let us conclude this subsection with an example such that Theorem 14.41 can be applied to get an 
SDP representation for Sv{f) while Theorem 14.21 can not. 

Example 4.6. Consider the convex set 5 = {(xi, X2) G : xi(x^ 4 X2) — xf — x^X2 — x| > 0}. The origin 
is a singular point on the boundary dS which is a quartic bean curve. The picture of this convex set is 
the shaded area bounded by the thick bean curve in Figure [31 After the perspective transformation, we get 

S = |(xi,X2) e : Xl > /(X2) := X2 4 i^j^j ■ /(X2) does not have structures required by Theorem 14.21 

Obviously T>' = (— cx),oo). We can check that /(X2) is convex over (— oo,cx)), so its epigraph epi(/) = 



n 



'qmod°°'' (xi — /) whicli Can be represented as 

(xi,X2) : 3y2,zo,zi,Xi > j/2 4 zo. 



Zo zi 1 - Zo 

Zl 1 - Zo X2 - Zl 

1 - Zo X2 - Zl 2/2 - 1 4 Zo 



>- 
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Figure 3: The singular convex set defined by xi{xi + x\) 



2 2 
X1X2 



> 



A direct SDP representation of 5* can be obtained by applying the inverse perspective transformation 



(.Ti,X2) : 32/2,20,2^1, 1 > 2/2 +20, 



zo zi xi - za 

Zi Xi - Zq X2 - Zi 

Xl - Zq X2- Zi y2- Xi+ Zq 



>- 



In the above, 1/2 + zq can be placed by one parameter. The plot of the above coincides with the shaded area 
in Figure [3l 



4.3. General case 

For the general case, we have the following result by applying Theorem 13.41 

Theorem 4.7. Assume V and Sf,{f) are both convex and have nonempty interior, and dim({i;i = 0} n 
dSf,{f)) <n-l. 

(i) If f{x) and every gi{x) are q-module concave over T) with respect to ij^^, then Sf,{f) = TZ^modif) ■ 

(ii) If f{x) and every gi{x) are preordering concave over V with respect to x'^~^ , then Sf,{f) — Ti^reif)- 

After one perspective transformation p, the singular point in Sx>{f) is mapped to one point at infinity of 
Sj){f), i.e., Sjj{f) itself does not have a point which is the image p(0). And the mapping p is smooth when 
Xl > 0. At any point x G Sx>{f) with xi > 0, the mapping p will preserve the singularity or nonsingularity 
at X. In this sense, the perspective transformation p will remove one or more singular points. Of course, the 
new convex set S-j=,{f) might have singularity somewhere else. In this case, we can apply some coordinate 
transformation to shift one singular point to the origin and then apply the perspective transformation again. 
So a sequence of perspective transformations might be applied. If there are finitely many singular points on 
the boundary, a finite number of perspective transformations can be applied to remove all the singularities. 
However, this approach might not work if there arc infinitely many singular points, i.e., the singular locus is 
positively dimensional. For instance, the convex set 



2\3 



{a;eM-^ : {1 - {xi - ly - x^) 



>0}, 



has a singular locus of dimension one. In this case, a finite number of perspective transformations is usually 
not able to remove all the singularies. 
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5 Some discussions 



Wc conclude this paper with some discussions and open questions. 



More general convex sets It is very natural to consider general convex sets of the form 



fk{x) > 0}, 



where fi{x), . . . , fk{x) are given polynomials or rational functions concave over the convex domain 2? = {x £ 
M" : gi{x) > 0, • • ■ ^gmix) > 0} defined by polynomials gi[x), ■ ■ ■ ,gm{x)- Note that 



So it suffices to consider each individual Sx){fi) separately. 

One interesting but unaddressed case is that the defining polynomials f i are concave in some neighborhood 
of Svifi, • • • , fk) but neither q-module nor prepordering concave over the domain V. In this situation, is 
it always possible to find another domain D' D Sx>{fi, • • • , fk) such that the fi is q-module or prepordering 
concave over V with respect to some other (p', q')l Or is it always possible to find a different set of defining 
polynomials for T) = {x E R" : gi{x) > 0, • • • , gm'{x) > 0} such that fi is q-module or prepordering concave 
over I? using new defining polynomials with respect to some different (p, q)7 This is an interesting future 
research topic. 

The separability in Positivestellensatz The rational function f{x) is concave over 2? if and only if 



for some sos polynomials r](x,u),a^^p,{x,u). Here fden is the denominator of f{x) which is nonnegative over 
V. When rj{x,u) = rii{x)r/2{u) is separable, we can choose p{x) = r/i{x) fden{x) and q = ?72('u)/Jg„(it), and 
then get an SDP representation for Svif) by following the approach in Section [3l However, in general case, 
is it always possible to find a factor i]{x, u) that is separable? Or what conditions make the factor 77(2;, u) to 
be separable? This is an interesting future research topic. 

Resolution of singularities In algebraic geometry [5], a well known result is that any singular algebraic 
variety (over a ground field with characteristic zero) is birational to a nonsingular algebraic variety. But 
the convexity might not be preserved by this birational transformation. Given a convex semialgebraic set 
in M" with singular boundary, is it is birational to a convex semialgebraic set with nonsingular boundary? 
Or is every convex semialgebraic set in M" equal to the projection of some higher dimensional convex 
semialgebraic set with nonsingular boundary? To the best knowledge of the author, all such questions are 
open. An interesting future work is to discuss how to remove the singular locus of convex semilagebraic sets 
while preserving the convexity. 

Acknowledgement The author would like to thank Bill Helton for fruitful discussions. 
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